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Abstract 

We examine the impact of threshold resummation for the inclusive hadronic production cross section 
of gluino pairs at next-to-next-to-leading-logarithmic accuracy, compared to the exact next-to- 
leading-order cross section and the next-to-next-to-leading-order approximation. Here, we apply 
formulas derived recently in the classical Mellin-space formalism. Moreover, we give the analytic 
input for the alternative momentum-space formalism and discuss the crucial points of the numeric 
implementation. We find that soft resummation keeps the hadronic cross section close to the fixed 
next-to- leading-order result. 



1 Introduction 



Within the search for new physics at the LHC experiments, one hopes to find evidences for heavy 
color-charged particles. These appear in many scenarios for physics beyond the Standard Model 
(SM). A prominent candidate is the gluino, which is a Majorana fermion and the superpartner of 
the SM gluon within various realizations of Supersymmetry (SUSY), as one of the most favorite 
models. Unfortunately, the direct search at ATLAS and CMS have only produced exclusion limit^ 
so far which, of course, depend on the model under consideration, and the assumptions on the SUSY 
particle spectrum. 

In order to separate possible signals from the SM background, a crucial theoretical input quantity 
is the inclusive hadronic production cross section (jpp^ggx- As known from standard Quantum 
Chromodynamics (QCD), the latter is given by a sum over the various partonic production channels 
where the partonic cross sections dij^gg are convoluted with their respective parton luminosity 
functions Lij. Explicitly, one has 

(^pp^ggx{s,m~,m~,ti'j,f4) = X] / dTLij{T,fij)aij^gg{TS,m~,m~,fxj,f4), (1) 

i,j=g,Q,g •^^™|/^ 

where the parton luminosities are themselves given by a convolution of the parton distribution 
functions (PDFs) 

Lij{T,fij) = dxi dX2 5{xiX2-T)fi/p{fi'j,Xi) fj/p{flj,X2) , (2) 

Jo Jo 

and we have introduced the hadronic center-of-mass (cms) energy s, the partonic cms energy s = ts, 
the factorization scale fip and the renormalization scale ^u,,. The gluino mass is denoted by rrig, and 
we assume mass degeneracy among the squarks flavors, therefore using a single scale niq. 

At leading order (LO) in perturbation theory, gluino pair production is driven by the partonic 
sub processes of gluon fusion gg — )• gg and quark-antiquark annihilation qq — )• gg, q = d,u, s, c, b. 
At next-to-leading-order (NLO), also the gq channel opens |3|. Close to the threshold however, its 
contribution is suppressed compared to gluon fusion and quark-antiquark annihilation. The full 
NLO result has been implemented in the public program Prospino [2]. 

As very well known, the cross section develops powers of large logarithms ln(/3) near the pro- 
duction threshold, where the velocity 



13 = Jl-Aml/s = (3) 
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of the produced particle pair goes to zero. The so-called threshold logarithms spoil the validity of 
the perturbative expansion in the strong coupling constant a^. However, they can be resummed 
systematically to all orders in perturbation theory, where the quantity as ln(/3) counts as order one. 
In this context, threshold logarithms are sometimes called soft logarithms and one talks about soft 
resummation. There are two approaches to soft resummation, which we will shortly discuss below. 
We will refer to them as Mellin-space formalism [5||6] and momentum-space formalism [7|j8]. 

Another difficulty arises through the exchange of soft gluons between the final state particles. 
At fixed order in perturbation theory, this gives rise to so-called Coulomb terms proportional to 
powers of as/ 13. Obviously, with decreasing /3, the above ratio becomes 0{1) and should therefore 
be resummed as well. This can be done in the framework of non-relativistic QCD (NRQCD). A 
joined soft and Coulomb resummation has been worked out in the momentum-space formalism [oj. 
In the Mellin-space formalism however, one would need to calculate the Mellin transformations 



^See Refs. |ll 2| for example. 
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of the NRQCD expressions analytically (or at least semi-analytically as a function of the Mellin 
moments), a problem which has not been solved so far. Therefore, Coulomb terms are included at 
fixed order only. 

For any process, soft resummation up to so-called next-to-leading-logarithmic (NLL) accuracy 
requires the knowledge of the color-decomposed Born cross section with respect to the SU{3)c- 
color configuration of the produced (s)particle pair. At next-to-next-to-leading-logarithmic (NNLL) 
precision, one needs the color decomposition of the NLO cross section near the threshold, where 
higher powers of /3 are skipped. For gluino pair production, the state of the art in the momentum- 
space formalism is a combined soft and Coulomb resummation up to NLL accuracy [10] . Within 
the Mellin-space approach, the NLL results of Ref. (lll|12| have recently been extended to NNLL 



Ref. 14 
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in Ref. 
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. In the context of NLL resummation, finite-width effects recently have been discussed 
It has been found that for T/rUg « 5%, the effect is of the order of the ambiguities due 
to higher-order contributions in the resummation formula. In this article, we work in the zero- width 
approximation which is very well justified if the squark and neutralino masses are about the order 
of the gluino mass. If one expands the NNLL threshold resummation formula up to (modulo 
the factor in the Born cross section), one reproduces the threshold logarithms up to next-to- 
next-to-leading-order (NNLO). One further obtains constant terms which are generically different 
at NNLO for the two approaches. However, a proper matching eliminates these constants and one 
is left with the threshold approximation (NNLOth), where, after factorization of the Born cross 
section, all NNLO constants are set to zero. At NLO, the constants are kept as these refer to the 
exact fixed order result. 

If the gluino is heavy (of the order of 1 TeV) , the threshold-enhanced terms give the dominant 
contribution at given order in perturbation theory. Combining the full NLO result with the 0{ag) 
threshold enhanced terms to the NNLO approximation (NNLOapprox)) one finds an increase of the 
cross section of about 20% with respect to the fixed NLO result flsl. The question of interest 
is how these findings change in the presence of NNLL resummation. Although the resummation 
formula for soft logarithms has been derived analytical, it has not been implemented yet. In this 
context, a general question is how to deal with ambiguities related to the choice of scales or the 
treatment of the Landau pole in momentum integrals. A different source of uncertainty comes from 
the unknown squark masses. However, the hadronic process is mainly driven by gluon fusion, which 
has only a weak dependence on rriq. As a consequence, the impact of the squark masses on the cross 
section is negligible to good approximation. In the threshold limit, the dependence on rriq actually 
vanishes. The impact of higher order contributions is estimated by a variation of the factorization 
and renormalization scale. As discussed in Ref. 13 , the main source of uncertainty comes from the 
shape of the gluon PDFs at high momentum fraction x, and the value of as{Mz)- The possible 
large discrepancies are not covered by the individual PDF errors, and give the main uncertainty 
when setting a reliable exclusion limit. 

In this article, we study the impact of NNLL soft resummation in the Mellin-space approach. 
In the sections [2] and [sj we briefly review the methods for soft (and Coulomb) resummation in the 
Mellin- and momentum-space formalism. Concerning the latter, we derive the color-decomposed 
hard function at NLO, which is needed for NNLL resummation. We discuss the ingredients of the 
potential function, which collects all results from NRQCD derived from the NLO non-relativistic 
potential. In the sections [4] and [5| we discuss technical aspects of the numerical implementation of 
the resummation formulas in the two approaches, and explicitly apply the Mellin-space formalism. 
The inclusive hadronic cross section is discussed in section [6] and we give our conclusion in section [7j 
An explicit formula for the matching of the momentum-space resummation formula onto the NNLO 
approximation is provided in the appendix. 



2 



2 Threshold resummation in MelUn space 

The traditional approach to threshold resummation has been invented in Ref . [5]|6] , see also Ref. (T^ 
[22] for further discussions. Resummation is performed in Mellin space after introducing moments 
with respect to the variable p = 4m~/s of momentum space. Neglecting the dependencies on fij 
and fij,, one has 
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ai}(m|) = M[<T.,(m|)](iV) = J dp p''^' a^jis,ml) , (4) 



where the threshold limit /3 — t- corresponds to — t- 00. The resummation of threshold logarithms 
is achieved by the formula 

= ^f/f9li9?^i{N + l)e^p[G,jM.N + ^)]+0{N-'ln^N). (5) 

All contributions from soft and collinear radiation exponentiate and are collected in the function 
Gij^i- To NNLL accuracy, it may conveniently be expanded according to 

G,j,i{N) = In(iV) • ^.(A) + 4. i(A) + a, g^iX) + ..., (6) 

where A^ = A^exp(7£;) with je denoting the Euler-Mascheroni constant, A = /3o InA^, and = 
as/(47r). The explicit expressions are calculated by a double integral over a set of anomalous- 
dimension functions. Some of these depend on the S'f7(3)c -color configuration of the final-state 
gluino pair, which we label by a capital index I. The color-summed partonic cross section for a 
given channel is simply given by 



Y^^ij,!- (7) 



The possible final-state color configurations are obtained by a decomposition of the initial color 
states into irreducible representations according to 

8x8 = 1, + 8, + 8a + 10 + 10 + 27,, (8) 
3x3 = 1, + 8, + 8a, (9) 

for gluon fusion and quark-antiquark annihilation, respectively. There are symmetric (s) and anti- 
symmetric (a) color states. The corresponding quadratic Casimir operators Ci, which show up in 
the color-decomposed partonic cross section, take the values 

Ci = {0,3,6,8} for I = {1, 8, 10, 27} (10) 

independent from the symmetry properties. The (color-decomposed) Born cross section ajj j factors 



out in Eq. (j5|), which is also true for the fixed order partonic cross section in the threshold limit o"*^ j. 
The matching constant g?- j collects all contributions independent of A^ and has to be determined 
order by order in perturbation theory by expanding the right-hand side of Eq. ([5]) in a, and 
matching onto o'fji- Coulomb corrections are not resummed by the above ansatz. They are formally 
treated as part of the hard scattering function, on which the exponent has to be matched. Therefore, 



one introduces a second matching term g^^j 23 as a power series in a,, which is actually a function 
of A'". If Coulomb corrections are not considered, g^y^ reduces to a factor one. 



^This produces all the threshold logarithms of the form ln''{N). 
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Note that starting from NNLO, one has non-relativistic corrections of non-Coulomb type, which 
cause the cross section to depend on the spin-configuration S of the produced heavy particle pair 
24-28 . These terms are also not treated by the exponent ^ and thus are fully hosted by the hard 



matching constant. As they also depend on the Mellin variable N, it is convenient to include them 
into g^^^-j^{N), although they do not arise from Coulomb exchange. Alternatively, one could also 
introduce a third matching constant g^j^'-"{N). The spin dependence is described by a parameter 
I'spinj which is zero for the gluino pair being in a spin singlet (S=0) configuration, and —2/3 for a 



spin triplet (S=l), see Ref. 24 . At the threshold, gluon fusion gives rise to a singlet, where quark- 
antiquark annihilation produces a triplet. All ingredients of the NNLL threshold resummation in 
Mellin space for gluino pair production can be found in Ref. |13|. 



3 Threshold resummation in momentum space 

An alternative approach, which features resummation in momentum space, has been invented in 
Ref. [7|[8] in the context of deep inelastic scattering. Here, soft and collinear radiation are described 
in the framework of soft-collinear effective theory. In Ref. [9j, this method has been adopted to 
heavy (s)particle pair production and extended to the summation of Coulomb corrections. The 
starting point is the observation that near the threshold the partonic production cross section 
can be factorized into a hard function H, a soft function W which collects soft fluctuations, and 
a potential function J which sums Coulomb exchange. Below the threshold, one has a discrete 
spectrum of bound-states. If these are taken into account, J should be expressed as a function of 
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the energy E = y/§ — 2mg relative to the threshold j9 

At fixed order in Og, the partonic cross section close to the threshold can be written as [9] 

aij^ggis,f,) = EE^§,i(^) / du;JiiE--)Wiiu;,f,), (11) 
s I -^^ 

where one sums over the spin and color configurations of the final-state particle pair. The spin 
dependence enters the hard and the potential function at NNLO first. In order to account for the 
different scales of the problem and to achieve resummation of threshold logarithms, one calculates 
the hard function at the hard scale and the soft function at soft scale fig- Then, one solves 
renormalization group equations for the hard and the soft function, which in the latter case sums 
the soft logarithms. The solutions are used to evolve H{fih) and W{fj,s) to the common scale fij, 
which is used for the convolution of the partonic cross sections with the parton luminosities. The 
general all order resummation formula reads [o] 

<7-|!^pg(s,/i/) = ^^H^jjiHh)Uij^i{mg,Hh,Hg,Hf) (12) 
s I 



where the summation of Sudakov-double logarithms is included in the evolution function C/jjj, 
which is given by 

/4m2\ "^"^^''"'^^^ fn'iy r 
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The function Sij^i{p,fis) is the Laplace transform of the MS-renormaHzed soft function Wi{uj,fi), 
which, to NLO accuracy, is given by 



{Ci + C, 



vr 
"6 



2Ci(p 



(14) 



12) is set to ?? = 2ar(/U^, p^) 



where we again define Og = as/(47r). The auxihary variable 77 in the Eq. 

after performing the derivative. It contains single logarithms [10] , whicE~can be seen by expanding 
rj = 4as(Ci + Cj) \n.{pg/ pj) + 0{a1). From the latter expression one also learns that 77 is negative 
and tends towards zero as Pg approaches pj. This behavior is preserved at higher orders in a^. 
As a consequence, the integration kernel (a;/(2mg))^^/a; has to be understood in a distributional 
sense [s], as will be further discussed in Sec. [5] 

The quadratic Casimir operators Cj and Ci depend on the color configuration of the initial and 
final states, respectively. For gluino pair production via gluon fusion we have Cj = Cj = Ca = 3, 
for quark-antiquark annihilation we have Cj = Cj = Cp = 4/3. With p^ < pj it hence follows 
that rj < 0. The function ar{i^, p) and the Sudakov exponent S{v^p) are given to NNLO precision 
in Eq. (86) and Eq. (87) of Ref. [30]. For NNLL resummation, the respective NLO expressions 
are sufficient. For i = j, these are functions of expansion coefficients F^ of the cusp anomalous 
dimension 



cusp 



En 

fc=0 



fc+1 



— Ci Tcusp ) 



(15) 



and coefficients of the QCD beta- function (5^- The latter are listed in Eq. (85) of the above article, 
whereas the cusp anomalous dimension is given for the quark-antiquark channel in Eq. (81). For 
gluon fusion, one just has to multiply these results by Ca/Cf, according to the Casimir scaling 



indicated in Eq. (15) above. For i 7^ j, one would replace 



cusp 



F'J = - (F 

^ cusp 2 ^ ' 



cusp ' cusp/ 



fc=0 



'cusp ' 



fc+1 



(16) 



Note also that the coefficients F^ coincide with the coefficients A, 



(0 



with I = A: -|- 1, used in the 
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Mellin-space formalism (see Ref. 
analogy to ar, by a replacement of 



31 



The functions afj j and d. 



'cusp with 'jY 



in Eq. (13) are obtained in 
+ 7j + IH, I in the first case, and with 



Hj,I - % 

(jf -|-7^)/2 in the second. The single particle anomalous dimension functions ■ji are related to soft 
radiation from the (mass-less) initial-state particles and are collected in appendix A of Ref. [32]. 
The anomalous dimension jh,i refers to soft radiation connected to a massive final-state particle 
in the color representation I. NLO results are given in Eqs. (3.31,32) of Ref. |33 . Here, one can 



also find expansion coefficients of 'jf, which are related to the evolution of the parton distribution 
function. Explicit expressions up to NLO for color-triplets and octets are given in Eqs. (D. 17-20). 

When taking results from the literature, care has to be taken with respect to the meaning of the 
parameter nj. Typically, nj denotes the number of active flavors, which contribute to the evolution 
of ag. So it does in the articles listed above, where the various coefficients of the QCD beta function 
and anomalous dimensions are collected. However, in some papers about gluino pair production at 
the thresholcQ nj denotes the total amount of quark (and squark) flavors, which, depending on the 
renormalization scheme, may differ from the number of active flavors. 



As discussed in Ref. |9j, the all-order resummation formula (12) does not depend on the hard and 
soft scale p^^ and Pg. However, as all ingredients are only known to finite order in a^, a truncation 



See Ref. |l3l 15| for instance. 
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of the perturbative expansion induces a residual dependence on these scales. A good perturbative 
behavior is obtained by a suitable choice for fih and fig- For the case of gluino pair production, 
these are fih ~ khm-g and fig = ksmg/S"^, with and kg being numbers of 0(1). The choice of a 
running soft scale ~ rrig/S'^ is also necessary to reproduce the threshold logarithms at NNLO (fixed 
order) [29]. In practice, the above choice becomes problematic when f3 tends zero. Moreover, as{fig) 
approaches the Landau pole. Therefore, one should introduce a parameter /3cut [29]. For /3 > /Jcut, 
a running soft scale is used, such that the perturbative expansion is not spoiled by large logarithms. 
On the other hand, for j3 < (3cut, the scale is fixed and kept in the perturbative regime. Note 
that a special treatment of the Landau pole is also required in the Mellin-space approach when 



performing the Mellin inversion (see Refs. [19,20 for instance). In the momentum-space approach a 
proper choice of /3cut should keep the ambiguities, which arise when matching the resummed to the 
fixed order cross section, as small as possible. For gluino pair production, this has been analyzed 
in Ref. [10| . Finally, we note that the choice of the soft scale actually determines what is being 
resummed. A dedicated discussion can be found in Ref. [34] . Within the Mellin-space formalism on 
the other hand, one has the implicit scale choices fig = rrig/N and fif^ = fXp and the two approaches 
are formally identical up to C?(l/iV)-terms (see for instance the discussion in Ref. [30|). 

Coming back to the momentum-space approach, one requires the potential function jf{E), 
which describes the exchange of soft gluons between the final-state particles, as well as higher order 
corrections coupled to these diagrams. At fixed order in perturbation theory, pure Coulomb correc- 
tions (which correspond to ladder diagrams) give rise to terms proportional to (os//?)". As already 
mentioned, they can be treated by methods of NRQCD using the LO non-relativistic potential, and 
have been found to be summed by a Sommerfeld factor A*-^ 35 . Starting from NLO, one also has to 
deal with terms multiplied by ag{ag/ P)^ . Their summation requires the non-relativistic potential 
at NLO |24l|26|. The LO and NLO Greens functions of the non-relativistic Schrodinger equation, 

9 , 29 , and the contributions to the potential function are 



G^°\{E) and G\^^-^{E), are given in Ref. 
obtained by taking two times the imaginary parts. For positive values of the energy E relative to 
the threshold, one finds the simple leading expression 



(1) 



4'\e) 





^ a) e{E) 



(17) 



ag 2tt^Di 




»(E), 



where the Sommerfeld factor is given by A*^(x) 
coefficient Di of the QCD potential 



x/ (exp(2;) — 1) , and we have introduced the color 



Di 



Ci-Ca = {-3,-3/2, 0, 1} for I = {1, 8, 10, 27} 



(18) 



Note that close to the threshold, where E <C 2mq, one has ^Jm-gjE = + 0{^) with /3 < 1. 



This gives rise to the expression A*-" (irag Di/ P) also described in Ref. 11, 13 , which makes the 



summation of l//3-terms evident. For Di < the potential is attractive, for Di > it is repulsive. 

The NLO Coulomb function ji^\E) can not be brought to a simple form. However, it can 
be expressed as a linear combination of dilogarithms and nested harmonic sums with a complex 
argument X{E) = ag{-Di) / {2^ -E /rrig) 



include ji^^ at fixed order in a,,. It reads 



29 



As a first approximation, it may be sufficient to 



4\E) = '^al2^'D, 



/3o In 



E 



+ 21n 



2mg 



ai 



(19) 
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with ai = 31/9 Ca — 10/9 ni . The choice of the scale is discussed in detail in Ref. [29] . The 
complete fixed order non-relativistic corrections at NNLO are obtained by adding to Eg. (|19[) the 
convolution of the C'(as)-Sommerfeld term with the NLO threshold approximation (without the 
Coulomb term), as well as a non-Coulomb contribution discussed below. The full result has been 
given in Ref. (24j first. 

If the potential is attractive, discrete bound states may develop below the threshold. Ignoring 
the gluino widths, these are described by 



n=l 



rrigasi-Di) 
2^ 



a. 



1 + ) 0{-D,) , 



(^ < 0) , (20) 



for S-wave production [9 ,29, where the bound-state energies are given by 

En - 



(21) 



The correction terms 5ri and ei stem from the NLO Coulomb Greens function |36] and can be 
found in Ref. [29], Eqs. (B.4) and (B.5). We stress that bound-state production (2 — )• 1 process) 
has different kinematics opposed to pair production (2—7-2 process) and therefore requires different 
fixed order formulas already at the Born level. The NLO cross section for gluinonium production 
has been studied in Ref. 
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As mentioned above, the non-relativistic potential also exhibits spin-dependent terms of non- 
Coulomb type starting from NNLO. Due to explicit results for the threshold enhanced part of the 
NNLO cross section fls], we know that these terms play a minor role for gluino pair production. 
They comply with around 5% of the total NNLO contribution to the partonic cross section. Thus, it 
is seems to be sufficient to include these terms at fixed order, which has also been done for top-quark 
pair production 29 . We follow the latter reference and define the factor 



A^c,i(^) 



Di{CA-2Di{l + v,^,^))e{E) 



(22) 



where the parameter fspm has been defined above. In summary, the potential function required for 
NNLL resummation is given by 



JhE) = 4'\e) (l + a2ln(/3)Agc,i(^)) +^r(^), 



(23) 



where In(^) k, mr,]. 

As last ingredient of the NNLL resummation formula, the hard function Hf- i{iJif^ is needed at 
NLO. At LO, it is related to the partonic Born cross section [9,10,29 




(24) 



where the approximations are valid in the limit s — t- 4m~. Near the threshold, the LO cross section 
is factored out from higher order contributions, where, apart from an overall factor (3 due to the 
Born term, all positive powers of /? are set to zero. Writing 



(1) 



r(2) 



(25) 



the required NLO hard functions are obtained by comparing the NLO fixed order threshold results 



with the NLO expansion of Eq. (12). For this purpose, it is sufficient to set = Ji^\ and we 
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require the various anomalous dimension functions at LO only. Using 7cusp = 4, = —2Ci, 

7^^*^^ = — 7^*^^ and E = rngP"^, we find in agreement with Eq. (D.3) from Ref. [oj 



~ res, exp / ^ 



2 



l+aJ-^^+4(Q + C,)ln2^=«2, 



(26) 



4(Ci + + C,)(4 + Lfg)) ln(8/32) + {d + Cj)L% 
2(Ci + (a + C,)(4 - L^/.))Lj^ + 12C7i + id + Cj) (^32 - 

(2Ci - m + Cj) ln(2) - (7^ + 7f )) L/h + (C. + C.^L^, + O (a^) 



where we have defined L/ 



./ft = ln(/i^V/x^-j, ana i.fg , 
scale has canceled between the soft function and the expansion of the evolution function. The terms 
proportional to Lf^ run the hard function i/^^^^ down to the factorization scale /ij, see also the 

discussion in Ref. [o]. 

The explicit evaluation of the soft and the potential function at NLO has produced all threshold 
enhanced terms proportional to ln^(8/3^) {k = 1,2) and which are easily identified within 
the fixed order cross section in the threshold limit, given in the Refs. [s] and 13 . However, the 
renormalization scale has been set equal to the factorization scale in the latter articles. Adopting this 
choice simplifies the right-hand side of Eq. (26) due to Lj/j = 0. From now on, we set fi^ 



and L/ 



The dependence on the soft 



as, and switch to the notation of Ret. 13 , where Lfg = L 



The NLO hard functions for 



gluino pair production in the momentum-space formalism then read 



HflM = 4C5-12Ci-C^(64-i^) 



2Ci-^CaH'^)]L,,-2CaLI 



H. 



(1) 



4Cf, 



12 Ca-Cf{ 64 



llvr^ 



(27) 



(28) 



2Ca-8Cf ln(2) - 2/3o + 6 C7j,]L^ - 2 CpLf^ 



in the MS scheme with = 5 active light (mass-less) quark fiavors. In the gluon fusion channel, only 
color-symmetric parts contribute near the threshold. Therefore, one has to sum over 1=1,8, 27. 
On the other hand, in the quark-antiquark annihilation channel, one only has contributions from 
the anti-symmetric octet to first approximation. The one-loop matching constants C^^j and Cf^g^ 
are given in Eqs. (39) and (40) of Ref. 



squared squark and gluino masses. 



13 



They are functions of the ratio 



m-fm- of the 



If one expands the resummation formula (12) up to 0{aj), one obtains some constants besides 



(2) 

the threshold logarithms at NNLO. In order to obtain the full NNLO hard coefficients H^ -^, a 
dedicated two-loop calculation is necessary to determine all constant terms. As these are not 
available at the moment, it is common to keep only threshold enhanced terms. Thus, if one wishes 
to match the NNLL resummation onto the NNLO approximation, one requires a cancellation of 
all constants in the resummation formulas at O(a^). Within the Mellin-space approach, this is 
achieved by a proper choice of g^j j in Eq. (j6]), see Ref. 13 . For an analogous treatment in the 
momentum-space approach, we need to compute the next order in Eq. (26). Keeping only the 
0{a^) contribution of the potential function, m|/3/(27r), and the hard function up to 0{as), the 
convolution with the NLO soft function and multiplication with the evolution function reproduces 
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the soft logarithms ln'^(/3) (n = 1,..,4) and their coefficients at NNLO if, and only if, we choose 
fig = ksTUgf^'^. The explicit fixed-order expressions can be found in Ref. 



13 



24 



The NNLO 

coefficient of the hard function has to be chosen such that the constants of the above expansion 
cancel at O(a^). Here, one should keep in mind that also the yet unknown NNLO soft function, if 
taken into account, produces constant pieces at 0{a'^). In our case, we just arrange the NNLO hard 
function is such a way, that all NNLO constant pieces cancel. This ensures a proper matching on the 
approximated NNLO cross section under usage of the NLO soft function. The explicit result is given 
in App. |A) The relative impact on the total hard function is rather small, typically at the percent 

(2) 

level. Finally, we mention that for NNLL+NNLOapprox accuracy, one can also set j = 0, if the 
matching procedure is performed at the NNLO level (see also the discussion in the next section). 
This differs from the choice in Eq. ( A 1 ) only in higher order terms beyond NNLL accuracy, and 



is applied in Ref. [29]. For NNLL resummation, all input functions to the resummation formula 
are needed at 0{as). Note that the NLO expression of the Sudakov exponent S{v,p) requires the 
coefficients of the cusp anomalous dimension and the QCD beta function up to NNLO (see Ref. [30] ) . 
All other anomalous dimension functions are required at NLO only. 



4 Implementation of the Mellin-space formalism 



In this section, we discuss the numerical implementation of the resummation formula ([sj). The 
following procedure is required in general [37] : First, as contributions of higher powers in /3 become 
numerically important when calculating the inclusive cross section, one has to match the resummed 
cross section to the full fixed order result. This is achieved by adding the latter, while subtracting 
all terms up to 0{as) (modulo the common prefactor a1) of the expanded NNLL expression. In 
order to perform the matching on the NNLOapprox partonic cross section, one needs to subtract 
all terms up to 0{a1) of expanded resummation formula (denoted by NNLL(2)) and add back the 



approximated NNLO result 29 



If one chooses to resum in Mellin space, an inverse Mellin transformation is required, which has 
to be done numerically. Due to omitted 1/iV contributions in the resummation formula ([s]), the 
matching onto the NNLO approximation will give a slightly different result compared to the NLO 
matching^ The resummed partonic cross section is given by 



^res,NNLO;, 



~NNLL,N .NNLL 2, N 
fj. . , — ! — CJ. . , — 



.NNLOapp rox 

'^ij^aa 



(29) 



The numerical Mellin inversion requires an analytic continuation of the resummed cross section 
to the complex plane. In our case, one simply assumes complex arguments of the logarithms. A 
more general discussion on that topic is provided in |38j[39j , and a good routine is included in the 
program ANCONT 38 . Here, the Mellin inversion is obtained by 



f{p) = - / dzIm[e^*p-^(^)M[/](c(z))] 
^ Jo 



(30) 



where Mellin- is identified with c{z) — 1, and c{z) = cq + ze . As the contour integral around 
the singularity at = is symmetric with respect to the x-axis, only the upper half is evaluated 
and multiplied by two. The path is chosen to be a truncated line with an angle $ close to tt. The 
parameter integral over z is divided logarithmically into 20 pieces, where each segment is performed 
by the 32-point Gauss formula. We choose the starting point cq = 1.4 to the right of the singularity 
at c{z) = 1, but left to the Landau pole singularity, which is excluded from the integration contour 



according to the minimal prescription presented in Ref. 19 , 20 



*Here, one has to emphasize that the numerically inverted NNLL(2) cross section does not exactly correspond to 
threshold terms NNLOth in momentum space, but comes close towards the threshold. 
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Figure 1: Partonic cross section for gluon fusion versus the energy above the production threshold normahzed 
to the gluino pair mass. The renormahzation and factorization scales have been set to fj, = m-g = 800 GeV. 
We plot the NNLO approximation which is exact up to NLO, and the NNLO threshold limit NNLOth which 
contains only threshold-enhanced contributions and NLO constants. Moreover, we show the NNLL resummed 
cross section and its expansion in up to second order, NNLL(2), as well as the resummed partonic cross 
section matched onto the NNLO approximation. First panel: Coulomb corrections are included at fixed order 
in the resummation formula. Second panel: Coulomb corrections are neglected during resummation but kept 
in the NNLO approximation. 



10 




Figure 2: Comparison of the resummed partonic cross section to the exact NLO result and the NNLO 
approximation. The scale is varied within the interval fi e [mg/2,2mg] around fi — rrig = 800 GeV. 



In Fig. [T| we plot the resummed partonic cross section in the gluon fusion channel against the 
energy above the threshold normalized to the gluino pair mass E/{2mg) = ^/S/{2mg) — 1 (solid 
black). The Coulomb corrections are included at fixed NNLO, and we set /i = m-g. We also plot 
the NNLO approximation (solid blue) and the threshold approximation NNLOth (dotted blue). 
Concerning the latter, all analytical input is given in Ref. |13|. As a check of the numerical Mellin 
inversion, we further show the inverted expansion of the soft resummation NNLL(2), which comes 
close to NNLOth for /3 < 0.5. The inverted soft resummation NNLL merges with its 0{a'^) expansion 
away from the thresholcQ For low velocities (/3 <C 0.1), the numerical quality of the Mellin inversion 
decreases. As this systematic error affects both NNLL and NNLL(2), one may expect a cancellation 
in the difference to some extent. One the other hand, /3 < 0.1 corresponds to an energy region 
E < W GeV for m-g ~ 1 TeV, which is of minor importance for the evaluation of the total hadronic 
cross section. In the second panel, the resummation is done by omitting all Coulomb corrections 
(these are still included in the NNLO approximation). In Fig. [2| we compare the resummed cross 
section to the exact NLO result, where we also plot the theoretical error related to scale variation in 
the interval /U G [mg/2, 2mg\. As evident from the figures, the NNLL soft resummation significantly 
depletes the partonic cross section compared to the NNLO approximation, and brings it down to 
the level of the NLO calculation for /3 > 0.1. Close to the threshold, the resummed cross section is 
enhanced compared to fixed NLO, but still well below the NNLO approximation. In the kinematical 
region /3 G [0.6,0.7], the resummed cross section, the NNLO approximation, and the threshold limit 
NNLOth cross each other close by. For higher velocities (/3 > 0.7), threshold approximations become 



untrustworthy. In our analysis we set both terms in the square brackets of Eq. (29) to zero, if (5 
exceeds the value 0.64. This choice gives rise to a smooth transition between the fixed-order and 
the resummed partonic cross section. Note that for the gluon fusion channel, the results are rather 
stable against a variation of the squark masses, which enter at NLO. In constructing Fig. [T| we 
have assumed as example m-g = 800 GeV and a common squark mass of 640 GeV. Choosing squark 
masses above 1 TeV will not change the conclusions drawn above. This is not true for the quark- 
antiquark annihilation channel. However, the hadronic production cross section is mainly driven by 
gluon fusion (about 99%). Thus, we will not have a closer look on gg-annihilation. 



^At high velocities (/3 > 0.95), the two results drift apart again. However, in this region the factorization ansatz 
in Eq. ([5| is not vahd and the resummation formula should not be used. 
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5 Implementation of the momentum-space formalism 



Within the momentum-space formalism, no inversion is required and one can remove the operator 
in Eq. ( 29 ) . However, the w-integral over the potential function in Eq. ( 12 ) requires an analytic 
continuation to negative values of t], where the integrand has to be understood in a distributional 
sense such that the integral is convergent. For NNLL resummation of gluino pairs we need to 
integrate 

1 



where 

J(0)( 



/(I) (a; 



jf\E-uj/2) ( 1 + 4a, 



oj V 2mg 



2*? e"2^-E 'n 



r(2r?) 



Ci[l + V'(2ry) + In 



2m~g 



vr 
24 



-\Ci + 2{Ci + Cj) ^{2ri)+\xi 
4^\E-oj/2) . 



'4>{2r]f - i)'{2ri) + 2 V)(2r/) In 



2mg 



2mg 



2m g 



In 



2m g 



+ {a + Cj) In^ 



2m g 



(31) 



(32) 



(33) 
(34) 



Here, we have defined ip^x) = ^e+^{x), and ■0 denotes the digamma function. Above the threshold, 
the integral runs from to 2E. As outlined in Ref. (8|[^ for instance, analytic continuation to 
negative values of r/ can be achieved by a replacement of the integration kernel by a so-called star 
distribution, which contains subtraction terms that render the integral finite for r] > —1. Explicitly, 
one applies 



2E 



dujf{E 



1 



UJ 



u \ 2mg 



2v 



2E 



du} 



/(i?-|)-/(s) + |m' 



UJ 

2mg 



2r? 



2E 



dujf{E--) 
'■^^duj 



ln( 



U! 



\2mg 



2r,- 



UJ 



fiE-^)-f{E) + '^f\E) 



In 



UJ 



+ 



m 

2r} 



Le 



1 

2ri 



r{E)E f 
2r]+l V 



Le 



27] +1 



+ 



f{E) f'{E)E 



2r] 



2mgJ \2mg 



UJ 



E 



m~g 



2ri 



2rj 



27] + I 



E 



2ri 



(35) 



(36) 



2E 



dujf{E-^) 



UJ 



UJ 

2m g 



2v 



2E 



+ 



duj 

UJ 

m 

27] 



'fiE-^)-f{E) + '^f'iE) 



In' 



UJ 



Li 



1 ^ 1 



V 



f'{E)E 

27]+ I 



UJ 



2mn/ \ 2m7T 



Li 



27] + I 



2r? 



Le + 



(2r/ + l)2 



(37) 



E 



2»? 
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with Le = ln{E/mg) and f{E — oj/2) being a smooth test function on the interval oj G [0, 2E]. For 
positive r], one can drop the star-brackets and the above relations become simple identities. For 
negative ??, diverging boundary terms at w = are removed by the star prescription. In deriving 



the fixed order expansion (26), we have assumed positive rj in order to perform the convolution 
over the LO potential function. We want to stress here that there are different methods for the 
analytic continuation. An alternative treatment, based on integration by parts, has been proposed 
in Ref. [29j. However, due to the identity theorem for holomorphic functions, the continuation to 
rj < is unique, as the prescriptions agree for positive rj^ 

A different ansatz for analytic continuation is required for the non- Coulomb corrections. Here, 
one may apply 



2E 



dL0jf\E 



2E 



duo 



^ 
(0) 



+ 



f{>{E 



a;. /'E-uj/2 
— m ' 
2^ 



rUg 



1 



CO 



oj V 2mg 



2ri 



(38) 



2' 



J<°'(E) + |J<»"(E) 



2mgJ 



V2 



mg 



2r] 



{Le -IE- ^{l + 2r?)) 



jf^'{E)E 
2r] + l 



{LE-lE-^{2 + 2r,)) 



m~gj 



for a 



where the terms in the last two lines have been evaluated for positive rj. See also Ref. [29 
slightly different treatment. The latter relations can also be applied to the logarithm in Eq. (19), if 
one includes the NLO corrections from the Coulomb potential at fixed order. 
The poles at rj 



-1/2 and rj 



-1 in Eqs. (35 38) are canceled by the overall factor l/r(2?7) 
in formula (31), see also the discussion in Ref. [8^^- The pole at r/ = is only canceled if there is 



no logarithm ln(ct;) contained in the integral. Otherwise, the prefactor l/r(2r/) is not sufficient to 
give a final result for — )• 0, as can be seen on the right-hand sides of Eqs. ([36]) and (37). However, 
this does not cause any trouble. The limit ry — )• corresponds to /i^ — )• //j, which, using a running 
soft scale as stated above, is achieved for large velocities (/3 ~ 1). However, the factorization ansatz 



in Eq. (11) is not valid in that kinematical regiorQ For instance, it is easy to see that some NLO 
terms in the soft function dominate the LO term, when r/ approaches zero. Thus, one should not 
apply resummation for (3 close to one. This will exclude the problematic case ry — t- 0. 

0,1,2 



Another complication arises due to the terms proportional to f{E) L'^ E"^^ with k 



m 



the formulas (|35||37|), and j[^\e) Le E'^'^ in Eq. (|38|, which in most cases diverge for — )• 0. In 



our case, only the repulsive potential J^y and its products with powers of \n.{E) evaluate to zero 
at the threshold. For r/ < —1/2, the convolution with the parton luminosity is ill-defined, see also 
the equivalent discussion for top-quark pair production in Ref. 29 . However, as outlined in the 



latter reference, the continuation to negative values of r/ can be extended to the convolution over 
the hadronic cms energy, producing a finite hadronic cross section. 

At this point, we want to note that something seems to be different compared to the treatment 
in Mellin-space. Suppose that we neglect all kind of Coulomb corrections. It follows that the fixed 
order partonic cross section has to vanish near the threshold due to the overall factor /3 in the born 
cross section. In Mellin space, this behavior is not spoiled by soft resummation. In the momentum- 
space approach in the absence of Coulomb corrections, the potential function reduces to the zeroth 
order in Eq.(17), which is proportional to ^J^Jmg k, (3. If this is multiplied by E'^'^ , one finds a 
singular behavior for £^ — )• and rj < —1/4 at the partonic level. Note that such a phenomenon does 
not occur in the discussion about deep inelastic scattering in Ref. [s] , simply because the momentum 
integral over runs from zero to Q^, with Q being the hard scale of the problem. 



''Here, I want to thank Pietro Falgari for clarifying discussions. 

^The same statement is also true for the corresponding ansatz of the Mellin-space approach in Eq. |5|, of course. 
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Figure 3: Inclusive hadronic cross section versus the squared cms energy. The error bars refer only to scale 
variation. Left panel: Fixed order for rrig = 800 GeV and rrig = 1 TeV. Right panel: NNLO approximation 
and soft NNLL resummation matched onto the latter. 



In any case one should state that resummation defines the partonic cross section as a distribution 
only, irrespective if Coulomb terms are resummed or nol|^ Therefore, there is no meaningful way 
to compare with the purely soft resummation in the Mellin-space formalism at the partonic level, 
where the threshold behavior mimics the fixed order in the absence of Coulomb terms, but adhoc 
assumptions have been made in order to deal with the Landau singularity. 

We close this discussion by noting that threshold resummation in general suffers from ambiguities 
dealing with the question of how to treat the Landau pole on the one hand, or such related to the 
choice of scales and the matching procedure on the other. The explicit implementation of the 
momentum-space formalism is passed over to the experts on that approach. 



6 Hadronic cross section 



In this section, we examine the impact of the NNLL resummation onto the inclusive hadronic cross 
section for the examples nig = 800 GeV and nig = 1 TeV. The current exclusion bounds suggest 
a minimum gluino mass of about 1 TeV [l|[2]. However, we want to stress again that these results 
depend strongly on the chosen PDF set and the value of as{Mz)- The libraries CTEQ6.6 |40| and 
MSTW08 (41] used in Ref. I] lead to similar results of the hadronic production cross section. The 
more recent sets ABMll (42 for instance differ from the latter at large parton momentum fraction x, 
which is of special importance for heavy (s)particle production. As discussed in Ref. [l3] , the usage of 
ABMll with as{Mz) = 0.1134±0.0011 compared to MSTW08 featuring a^(Mz) = 0.1171±0.0014 
reduces the hadronic production by a factor of about two. Concerning the relative effect of the 
NNLL resummation compared to fixed NLO or NNLOapprox; the actual choice of the PDF set will 
only play a minor role. In our analysis, we always use ABMll NNLO, also for the LO and NLO 
cross section which we plot for illustration of the NNLO threshold effects. We further assume a 
fixed ratio nig /nig = 4/5 giving rise to the values nig = 640 GeV and nig = 800 GeV. As mentioned 
above, this choice has also no significant impact on the result. A variation of the squark masses 
would affect the inclusive hadronic cross section at the percent level only. 

*Here, I want to thank Martin Beneke and Christian Schwinn for clarifying discussions. 
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The theoretical error due to higher order corrections is estimated by a variation of the scale 
jJL = jJLf = within the range jj, E [mg/2, 2mg]. Within the momentum space approach, one further 
has the possibility for a variation of the soft scale in order to account for the ambiguities within 
the soft resummation, as well as an independent variation of the hard scale. In principle, the latter 
is also possible in the Mellin-space approach. For illustration of the perturbative accuracy, we will 
not include the PDF errors. Plots showing these can be found in Ref. fls]. The convolution of the 
(resummed) partonic cross sections with the parton luminosity functions according to Eq. ([T]) is 
performed by a VEGAS integration routine using LHAPDF grid files [43] for the PDF sets. 

In the left panel of Fig. |3| we plot the fixed order inclusive cross section as a function of the 
hadronic cms energy ^/s. For nig = 800 GeV and ^/s = 7 TeV, we observe the K factors -HTnlo = 
<7NLo/o-LO = 1-55 and -f^NNLO = o-NNLOapprox/o"NLO = 1-29, at 14 TeV we have ZCnlo = 1-50 and 
-f'^NNLO = 1-21. For nig = 1 TeV a similar picture emerges, with Knlo = 1-54 and i^NNLO = 1-33 
for ^/s = 7 TeV, and -STnlo = 1-48 and -f^NNLO = 1-23 at 14 TeV. In the right panel, we compare 
the NNLO approximation to the NNLL resummation, matched onto the latter. In both cases, the 
K factors Knnll = CTres.NNLOapprox/^NLO are slightly below one over the whole range of hadronic 
cms energies. This is a direct consequence of the partonic result shown in Fig. [2j where after 
convolution with the parton luminosities, the corrections due to resummation arise mainly from the 
region /3 G [0.1,0.6]. 

Concerning the variation of the factorization and the renormalization scale, we find a rather large 
error for the NNLL resummation, even exceeding the uncertainty of the fixed NLO computation. 
This can also be understood from the partonic cross section shown in Fig. [2] In the region of 
relevance, the NNLL uncertainty is slightly larger then its NLO counterpart. At the same time the 
central values are lower, hence giving rise to a stronger relative impact of the scale variation. We 
conclude that for the case of heavy gluino pairs, soft resummation suggests that the NLO description 
already gives a good theoretical estimate. An outstanding question is if this behavior is kept after 
Coulomb resummation, which is likely to be answered within the momentum-space approach soon. 



7 Conclusions 

In this article, we have discussed the effect of NNLL soft resummation on the inclusive cross section 
for gluino pair production at the LHC by the use of the Mellin-space resummation formalism. Here, 



we have used analytic results of Ref. 13 , and performed the Mellin inversion numerically by the 
use of standard methods. We have found that soft resummation compensates the enhancement of 
the NNLO soft logarithms compared to the fixed NLO result. This in turn suggests that the NLO 
computation already provides a good theoretical input for constructing exclusion limits at the LHC 
or, hopefully, to construct confidence regions for future evidences. At this point we want to stress 
again that the main source of uncertainty is related to the non-perturbative input, when different 
PDF sets and initial values for as{Mz) are applied. As discussed in Ref. |13| , the individual PDF 
errors are not sufficient to cover the resulting discrepancies. 

We further calculated the hard function required for a joined soft and Coulomb resummation in 
the momentum-space approach. Here, we recaptured the relevant points concerning the convolution 
over the hadronic cms energy, which are given in greater detail in the literature. Finally, we state 
that threshold resummation always suffers from ambiguities, no matter which approach is applied. 
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A Matching to the NNLO approximation 

In this appendix, we give the relevant part of the NNLO hard function, which may be used for the 
matching onto the NNLO approximation. The specific choice given here sets all NNLO constants 



to zero within the resummation formula (12). For the generic scale choice /i^ = ks'iTig/S'^, we obtain 



^(2)nratch ^ _ci^^^f^48Ci{l - ln(2)) + 2(Q + C,)(64 - 961n(2) + 721n2(2) - yvr^) j ([A}1) 

+1UC^{1 - ln(2)) + 2{Ci + Cj)Ci(384 - 4801n(2) + 2161n2(2) + Il7r2ln(2) - 2271^) 

+2{Ci + Cj f {512 - 7681n(2) + 5761n2(2) + Mtt^ ln(2) - SSvr^ - — vr^ + — vr^) 

+L^{-C;fi(8Ci + 16(Q + Cj){2 - 31n(2))) - 4C7i(33 - 2ni - Snij) (l - ln(2)) 

+24C72 - {Q + Cj)(^ - Jci(l3^2 ^ i44in(2)(l - 21n(2))) + ^(-272 + 4081n(2) 
V o 3 9 

-216 ln^{2) + llvr^) + Hij (-32 + 48 ln(2) - 36 ln'^{2) + y vr^) - 268 ln(2) - 528 ln(2) 
+3961n2(2) + 127T^ ln(2) - ^vr^) - ^(Ci + Cjf {768 - 9601n(2) + 10081n2(2) 



-4321n(2)3 + 6l7r2 ln(2) - 48^^ _ 336^(3)) + 7^ ^ + if'^ + tJ,?} + Ll[-ACl^{Ci + C, 

2 / 199 2 

-Ci(ll - -m - 2Kij) + 2Cl + (C, + Cj)[-— + 4Ci(3 - 21n(2)) + -n,(l7 - 181n(2)) 

+4Kij(2 -31n(2)) +661n(2) +7r2) + (Ci + Cj)2(-32 + 641n(2) - 241n2(2) + yvr^)} 

+L3 {(C, + C7,)(-^ + 2Ci + + K,,) - 4(Q + C,f ln(2)} + \{C, + Cjflf^ , 

ln(A:,){-(7f + 7^ - 7^(1^^ + 4Ci(33 - 2ni) (l - ln(2)) + 24^^(1 - ln(2)) 



+ 

+ 7,{Ci + Ci)( + ^Ci{672 - 1008 ln(2) + 6481n(2)2 - 357r2) - ^n/(272 - 4081n(2) 
2 V 3 3 9 



1.^ ^ ,/3184 2_ , , , , 2 

+216 ln(2)2 - llvr^) - 1592 ln(2) + 792 ln(2)2 + 247r2 ln(2) vr^ 

3 

+^{Ci + Cjf (1152 - 1728 ln(2) + 12961n(2)2 - 6481n(2)3 + 1057r2ln(2) - 707r2 - 336C3) } 

+ ln\ks)[2Ci + Ci(ll - + ^{d + C,)(^ + 8Ci(7 - 91n(2)) - n;(^ - 81n(2)) 

-1321n(2) - 27r2) + ^(Q + Cjf{576 - 8641n(2) + 6481n(2)2 - 35ir^)^ 

+ ln^{ks){^{a + C,) (y + 4Ci - ^n;) + 4(C, + C,f{2 - 31n(2))} + \n\ks)^{Ci + Cj)^ , 

where Hgg = and Kgg = —2(3o + 6Ci;'. The relative correction to the NLO hard function is about 
a few percent only. 
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